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We study the dynamics of an optical mode in a cavity with a movable mirror subject to quantum 
Brownian motion. We study the phase noise power spectrum of the output light, and we describe the 
mirror Brownian motion, which is responsible for the thermal noise contribution, using the quantum 
Langevin approach. We show that the standard quantum Langevin equations, supplemented with 
the appropriate non-Markovian correlation functions, provide an adequate description of Brownian 
motion. 



I. INTRODUCTION 



The mechanical interaction between a moving mirror and a radiation field has been an important topic for the study 
of very high precision optical interferometers in which radiation pressure effects cannot be ignored. This interaction 
is at the basis of the interferometric detection of gravitational waves, where the tiny displacement of a mirror can be 
detected as a phase shift of the interference fringes . Another interesting application is the atomic force microscope 
, where an image of a surface at atomic resolution is obtained from the measurement of the force between the 
surface and a probe tip mounted on a microcantilever. 

A cavity with a movable mirror is of interest also for cavity QED studies, which usually involves the quantum 
coherent interaction between high-Q cavity modes at low photon number and single atoms. In this case, the atomic 
degrees of freedom are replaced by the motional degree of freedom of the movable mirror. Interesting quantum effects, 
as the generation of sub-Poissonian light || , of Schrodinger cat states of both the cavity mode [Q and even of the 
mirror [|| have been already illustrated. 

In these applications one needs a very high resolution for position measurements and a good control of the various 
noise sources, because one has to detect the effect of a very weak force. As shown by the pioneering work of Braginsky 
H , even though all classical noise sources had been minimized, the detection of gravitational waves would be ultimately 
determined by quantum fluctuations and the Heisenberg uncertainty principle. Quantum noise in interferometers has 
two fundamental sources, the photon shot noise of the laser beam, prevailing at low laser intensity, and the fluctuations 
of the mirror position due to radiation pressure, which is proportional to the incident laser power. This radiation 
pressure noise is the so-called "back-action noise" arising from the fact that intensity fluctuations affect the momentum 
fluctuations of the mirror, which are then fed back into the position by the dynamics of the mirror. The two quantum 
noises are minimized at an optimal, intermediate, laser power, yielding the so-called standard quantum limit (SQL), 
which coincides with mean square fluctuations of the harmonic oscillator ground state AqsQL = \ffi/2mujs {uos is 
the mirror oscillation frequency). Real devices constructed up to now are still far from the standard quantum limit 
because quantum noise is much smaller than that of classical origin, which is essentially given by thermal noise. In 
fact, present interferometric gravitational wave detectors are limited by the Brownian motion of the suspended mirrors 
0, which can be decomposed into suspension and internal (i.e. of internal acoustic modes) thermal noise. Therefore it 
is very important to establish the experimental limitations determined by thermal noise and recent experiments [|]f| 
have obtained interesting results. With this respect it is also important to establish which is the most appropriate 
formal description of quantum Brownian motion. In fact, even though the classical understanding of the phenomenon 
is well established, relying on Langevin or Fokker-Planck equations fLO] , its quantum generalization is still the subject 
of an intense debate (see jy] and references therein). In particular, the recent paper by Jacobs et al. JlJ] has shown 
that the standard description of quantum Brownian motion, which is the straightforward generalization of the classical 
case fl3|| , gives an inadequate description since it generates a non-sensical term in the optical phase noise spectrum 
in the case of a cavity with a movable mirror. The authors of jllj adopt therefore a corrected quantum Langevin 
equation, based on the Diosi master equation and suggest that, even if it is quite challenging, the corresponding 
modifications of the phase noise spectrum could be revealed using miniature high-frequency mechanical oscillators 
and ultra-low temperatures. In the present paper we shall reconsider the same system, i.e. a driven cavity with a 
movable mirror, and shall show that the inadequacy shown in jl2| has to be traced back to the inadequacy of the 
quantum noise commutation relations and correlation functions which are dictated by the standard Brownian motion 
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master equation. We shall see that, differently from the master equation approach, a consistently applied quantum 
Langevin equation |H| provides a flexible approach, valid at any temperature and therefore also in the fully quantum 
regime of very low temperatures. This however does not mean that the quantum Langevin equation approach is 
generally superior than the master equation approach, but simply that in the case under study, which is a linearized, 
non-markovian problem, the quantum Langevin description is more convenient and powerful. 

The outline of the paper is as follows. In Sec. II the appropriate quantum Langevin equations for a Brownian 
particle are derived starting from the usual model based on the coupling with a reservoir of harmonic oscillator, and 
its consistency is shown. In Sec. Ill the quantum Langevin approach is applied to the case of a cavity mode with a 
movable mirror and the homodyne spectrum of the reflected light, showing the thermal and quantum fluctuations of 
the mirror, is studied. Sec. IV is for concluding remarks. 



II. THE DYNAMICS OF THE SYSTEM 



The system studied in the present paper consists of a coherently driven optical cavity with a moving mirror. This 
opto-mechanical system can represent one arm of an interferometer able to detect weak forces as those associated with 
gravitational waves Q| or an atomic force microscope j|] • The detection of very weak forces requires having quantum 
limited devices, whose sensitivity is ultimately determined by the quantum fluctuations. For this reason we shall 
describe the mirror as a quantum mechanical harmonic oscillator with mass m and frequency wg. The optomechanical 
coupling between the mirror and the cavity field is realized by the radiation pressure. The electromagnetic field 
exerts a force on the movable mirror which is proportional to the intensity of the field, which, at the same time, 
is phase-shifted by 2kq, where k is the wave vector and q is the mirror displacement from the equilibrium position. 
In the adiabatic limit in which the mirror frequency is much smaller than the cavity free spectral range c/2L (L is 
the cavity length) |l6f| , one can focus on one cavity mode only because photon scattering into other modes can be 
neglected, and one has the following Hamiltonian 

H = hw c tfb + — + -muj 2 s q 2 - h—qtfb + ihE - be iuJot ) , (2.1) 

2t7£ 2 L 

where b is the cavity mode annihilation operator with optical frequency oj c and E describes the coherent input field 
with frequency cl>o ~ lo c driving the cavity. The quantity E is related to the input laser power P by E = ^Pj c /huJo, 
where j c is the cavity decay constant due to the input coupling mirror. Since we shall focus on the quantum and 
thermal noise of the system, we shall neglect all the experimental sources of noise, i.e., we shall assume that the driving 
laser is stabilized in intensity and frequency. This means neglecting all the fluctuations of the complex parameter E. 
Including these s upp lementary noise sources is however quite straightforward and a detailed calculation of their effect 
is shown in Ref. . Moreover recent experiments have shown that classical laser noise ca n be made negligible in the 



relevant frequency range |||9|. The adiabatic regime u>s "C c/2L we have assumed in Eq. (2.1) implies u>s "C oj c , and 
therefore the generation of photons due to the Casimir effect, and also retardation and Doppler effects are completely 
negligible. 



The dynamics of the system is not only determined by the Hamiltonian interaction (2.1), but also by the dissipative 
interaction with external degrees of freedom. The cavity mode is damped due to the photon leakage through the 
mirrors which couple the cavity mode with the continuum of the outside electromagnetic modes. For simplicity we 
assume that the movable mirror has perfect reflectivity and that transmission takes place through the other "fixed" 
mirror only (see Fig. 1 for a schematic description of the system). The mechanical oscillator, which may represent 
not only the center-of-mass degree of freedom of the mirror, but also a torsional degree of freedom as in qJ, or an 
internal acoustic mode as in Q , undergoes Brownian motion caused by the uncontrolled coupling with other internal 
and external modes at the equilibrium temperature T. 

The dissipative dynamics of the optical cavity mode is well described by the so-called vacuum optical master 
equation 17 



P=Y (26p& t - tfbp - ptfb) , (2.2) 

for the time evolution of the density matrix of the whole system p. In fact, the mean thermal number of photons 
at the optical frequency lu c is e xtremely small and thermal excitation is therefore completely negligible. The time 



evolution generated by Eq. (2.2) presents no ambiguity. In fact it is of Lindblad form |]18| and therefore it preserves 
the positivity of the density matrix. Moreover it is completely equivalent to the time evolution for the operators in 
the Heisenberg representation driven by the following quantum Langevin equation 
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6(*) = -yK*) + V 5 Pin(*) J (2-3) 

where bi n (t) is the input noise operator associated with the vacuum fluctuations of the continuum of modes outside 
the cavity, having the following commutation relation 

~bin(t),bl(t')} =5{t-t') (2.4) 

and correlation functions 

<M<)MO) = (bt(t)b in (t')) = (2.5) 

(b in (t)bi(t'))=5(t-t'). (2.6) 

The description of the quantum Brownian motion of a massive particle in a potential is instead not so well estab- 
lished. The standard Brownian motion master equation (SBMME) has been first derived by Caldeira and Leggett 
|[13fl in the high temperature limit and reads 



/'(/'I = ~ [H S ,p] - ^ [q(t),{p(t),p(t)}] - ^ [q(t), [q(t),p(t)}] , (2.7) 



where Hg is the uncoupled particle Hamiltonian, p is its momentum and rjq is the friction force. This master equation 
is the direct generalization of the Fokker-Planck equation for the classical Brownian motion but, since it is not of 
Lindblad form, it has the drawback that it does not ensure the positivity of the density operator This fact 

has stimulated many authors who have amended the SBMME with additional terms so to cast it into the Lindblad 
form (ll] . These corrected master equations preserve the positivity of the density matrix but also them are 
valid in the high temperature limit only as the SBMME, because they necessarily provide a Markovian description 
of Brownian motion which instead becomes highly non-Mar kovian in the low temperature limit ]2^] . In fact, in this 
limit, the reservoir correlation time is no more negligible because it is essentially determined by the "thermal time" 
r T = H/kT. 

An alternative description of quantum Brownian motion is provided by the quantum Langevin equations for the 
Heisenberg operators q(t) and p(t), which, in analogy with the classical case, should read 

q{t)= l -\H s ,q(t)\ (2.8) 

p(t) = UH s ,p(t)]--p(t)+ttt), (2-9) 
a m 

where £(i) is a Hcrmitian noise operator with correlation function 

(£(*)£ (O) = 2vkT5(t - t') . (2.10) 

One would expect that these equations give correct results at least in the high temperature limit, where the classical 
limit should be recovered. Instead Ref. has shown that they give inconsistent results even in this limit because 
they do not preserve the commutation relation [q(t),p(t)] — ifi and they yield a spurious term in the phase fluctuation 
spectrum. In fact, in spectral measurements, the Fourier transform of correlation functions of the form G(t) = 
(R(t)R(t + t)) are measured, where R{t) is an appropriate output field. This correlation function depends only on r 
because of stationarity, and moreover it is an even function of r because R(t) commutes with itself at different times. 
This implies that the observed spectrum has to be an even function of the frequency u>, while the adoption of the 



quantum Langevin equations (|2.8j)-(2.£) yields a term which is an odd function of u> 



In Ref. jl2j these inadequacies of th e qu antum Langevin equations (2.8)-(BJj) have been traced back to the fact that 



they are equivalent to the SBMME (2/7), which is not of Lindblad form and therefore does not preserve positivity. 
For this reason they consider the amended master equation of the Lindblad form proposed by Diosi in Jl4| , and derive 
the set of quantum Langevin equation equivalent to it, 

q{t) = jr[Hs,q(t)] + <t) (2-H) 
p(t) = ^[H S)P (t))-^P(t)+m, (2-12) 
having the additional noise term e(t). The corresponding correlation functions are 



3 



(£(*)£ (*0) = 2rikT5{t - t') (2.13) 

(e(t) e (f)> = ^(t-0 (2.14) 

«(t)e(0> = -i^(*-0 (2-15) 

(e(t)e(t')>=*^(t-0- (2-16) 
2m 

It is then possible to see that the phase noise spectrum associated with this dynamical description of quantum 
Brownian motion has no spurious term and that it is an even function of the frequency, as it must be. However, the 
approach of Ref. can be questioned for two reasons. First of all the master equation of Ref. [Q has been derived 
by Diosi with heuristic arguments and in a later paper p0| Diosi himself corrected it by considering a more rigorous 
medium temperatures extension of the SBMME. The master equation of Ref. |^(J would lead to a different set of 
quantum Langevin equations; more generally speaking, the Lindblad form condition does not uniquely determine the 
master equation, and therefore the form of the quantum Langevin equations neither. Furthermore, the added noise 



term e(t) in Eq. (2.11) is not present in the classical case and it has an unclear physical origin. 

For this reason we reconsider here the problem, assuming a different starting point. Most of the derivations of 
the Brownian motion master equations are based on the independent oscillator model for the reservoir, whose quite 
general validity has been extensively discussed in p3| . Therefore, rather than first deriving the master equation from 
this reservoir model and then considering the quantum Langevin equations associated to it, we derive the quantum 
Langevin equations directly from the reservoir oscillator model, as it is shown in |TE , p4j ]. 

Let us neglect for the moment the presence of the cavity mode and consider a particle undergoing Brownian motion, 
with Hamiltonian H$. The reservoir is described by a collection of independent harmonic oscillators with frequency 
u)j, couplings A:,-, and whose canonical coordinates qj and Pj have been appropriately rescaled. The total system 
Hamiltonian is |l|,|l|,|4} 

H s + lj2^- k ^ 2+LU h")- (2-17) 

3 

The quantum Langevin equations can be obtained from the Heisenberg equations for q(t), p(t) and the reservoir 
annihilation operators Oj(t) = (uijqj(t) + ipj(t)) / W2hLJj, 

q {t)= l -[H s ,q{t)] (2.18) 
p(t) = l - [H s ,p{t)\ + ]T ki ( Pj (t) - k 3 q(t)) (2.19) 

3 

dj (t) = -iUjOj (t) - k s J^jrQ® • (2.20) 

If we integrate the equation for a,j(t) starting from the initial time to, we get 

I — pt 

a j (t) = e- iu ' ( - t - t ^ai(t )-k j J 1 ^ / Me^^-^qW), (2.21) 
which, using integration by parts and the fact that q = i [Hg, q] = p/m, can be rewritten as 

ai(t) = iki^m + e-^-o^o) - ^^(^) - ^-^'Mt'). (2-22) 

Then we replace ( 2.2 2| ) in the equation for p(t), ( 2.19 ), which becomes 



p(t) = j [H s ,p(t)] + Q(t) - V k) cos [ Uj (t - t Q )} q(t Q ) - f dt'Y] k) cos [ Uj (t - t')] (2.23) 



3 Jto 3 



where we have defined the reservoir operator 



4 



Q(t) 



,iwj(t-t ) t 



am) 



-iWj{t—to) 



i(to)) 



(2.24) 



As it is well known, the irreversible properties of the reservoir are obtained only when an infinite number of oscillators, 
distributed over a continuum of frequencies, is considered. The continuous limit has to be performed according to the 
following prescription 



dui k 2 (uj) 



dn 
du> 



2i) 



duj ■ ■ ■ 



(2.25) 



where dn/dio is the oscillators density, n is just the friction coefficient, and fl is the frequency cutoff of the reservoir 
oscillator spectrum. When the continuous limit is considered, the Heisenberg equations for the Brownian particle 
become 



m = & 

m 



p(t) = - [H s ,p(t)] + Q(t) - 2 v 5(t - t )q(t ) - 277 / dt'S(t - t') 1 - 

h Jto m 



where we have defined the following function 



8{t) 



dio cos(wi). 



(2.26) 
(2.27) 

(2.28) 



Eqs. ( [2.26 ) and (2.27) become identical to the usual Langevin equations (2.£)- fl2.9| ) when the usual assumption of 
a reservoir dynamics much faster than that of the Brownian particle is made. This means making a coarse-grained 
description in time equivalent to assuming the infinite cutoff limit Q — > 00, under which the function S(t) becomes 
a Dirac delta function. The reservoir operator Q(t) in Eq. ( |2.27 ) plays therefore the role of the random Langevin 
force, commuting with a generic system operator evaluated at the initial time to, for every value of t. This fact 
suggests to interpret Q(t) as the input noise of the system and the above quantum Langevin equations in terms of 
the input-output formalism developed by Gardiner and Collett J25|. However the interpretation of Q(t) as an input 
noise must be mad e w ith care, because its commu tatio n relations are different from those of t he ty pical input noise 
operators (see Eq. (2.4)). In fact, using definition ( 2.24 ) and the continuous limit prescription ( 2.25| ), one derives the 
following commutation relation 



Q(t),Q(t') 



2ihr)j t ~5(t - 0, 



(2.29) 



which is not a delta function in (t — t'), even in the coarse-grained time limit f2 — * 00. With this respect, it is 
interesting to consider also the correlation function of the noise operator Q(t), which is generally defined as a trace 
over the reservoir degres of freedom, 



(Q(t)Q(t')) B = tr B {Q(t)Q(t')Ro}, 
where Rq is the density operator of the reservoir at thermal equilibrium at temperature T, 



kT \ 



(2.30) 



(2.31) 



Using again Eqs. (|2.24j ) and ( [2.25 ), one gets 

{Q{t)Q{t'))B = ^-{Mt - + - *')}, 



with 



FrW) — I dio uj cos(wr) coth ( — — I 
Jo \2kTJ 



Mr) 



du> to sin(ojr) = 7r — S(r), 
dr 



(2.32) 

(2.33) 
(2.34) 
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The antisymmetric part, corresponding to Ti, is a direct consequence of the commutation relations ( [2.29 ) and, as we 
have seen, is never a Dirac delta, while the symmetric part, corresponding to J 7 ,., explicitely depends on temperature 
and becomes proportional to a Dirac delta function only i f the high temp e ratur e limit kT ^> Kl first, and the 
infinite frequency cutoff limit Q, — > oo later, are taken. Eqs. ( |2.29|) and (|2.32|) -( 2.34 ) show the non-Markovian nature 
of quantum Brownian motion, which becomes particularly evident in the low temperature limit p2[ . Therefore, 
assuming the independent oscillator model for the reservoir, which is the usual starting point for th e deri vatio n of the 
Brownian motion master equation, we have derived the exac t quantum Langevin equations ( 2.26| ) and (2.27), which 



reduce to the usual quantum Langevin equations (2.8)-(2.9) in the limit f2 — > oo. However, these equations must 



not be used as usual input-output equation s, as it is implicitely dictated by the S BMM E of Eq. (2.7) (see Ref. Jl2[) 



but the appropriate commutation relations ( 2.29 ) and correlation functions ( 2.32 )-( 2.34 ) must be used. Therefore the 
inadequacies found in Ref. |12| are not due to the form of the Langevin equations but only to the inappropriate form 
of the noise correlation function dictated by the SBMME. It is also important to stress that our Langevin equation 
description of quantum Brownian motion is more general than that associated with a master equation approach, 
because it is valid at all temperatures and it does not need any high temperature limit. 

Another criticism to the standard quantum Langevin equations presented in Ref. |l2T| is that they do not preserve the 
commutation relations of system operators. This can also b e trac ed ba ck to the inappropriate form of the commutation 
relations of the noise operator Q(t). Actually, Eqs. ( 2.26 ) and ( 2.27 ) are exact and, because of the unitarity of the 



time evolution, they mantain the initial comm utati o n ru les for the system operators. This pro perty i s pres erved also 
by the standard quantum Langevin equations ( |2.8| )-(f2~S|), which are the £1 — > oo limit of Eqs. fl2.26 )-( p^27 ), provided 



that the correct noise commutation relation is used. Let u s see this fact in detail. If we do not restrict to the usual 
condition t > t , under the fl — > oo limit, Eqs. ( [2.26| ) and (2.27) have to be written as (see Appendix) 



q(t) = p(t)/m 



[H S ,P] + Q(t) - 2r]S{t - t )q(t ) - ^-S(t - t )p(t) , 



(2.35) 
(2.36) 



where S(t) is the sign function (defined so that S(0) = 0). 

If we consider for example the commutator between q(i) and p(t) , differentiate it with respect to t and use Eqs. ( 2.36 ) 
we obtain 



-[q(t),p(t)]=-+S(t-t )[q(t),p(t)] + [q(t),Q(t)] . 



(2.37) 



In order to solve this equation we need the commutator between q(t) and Q(t). More in general we consider both 
quantities X(t) = [q(t),Q(t')] and y{t) = [p(t), Q(t')] as functions of t, with t' an independent parameter. If for 
simplicity we restrict to the case of interest here of a harmonically bound Brownian particle H$ — m,L0 2 s q 2 /2, it is 



possible to obtain the equation for X{t) and y(t) using Eqs. (2.36) 

X(t) = y(t)/m 

y{t) = -muj%X{t) + 2ihr)5'(t - t') - ^-S(t - t )y(t) , 



(2.38) 
(2.39) 



where w e hav e used Eq. ( 2.29 ) in the fl — ► oo limit and the initial condition X(to) = 0. Using Eq. ( 2.36| ), the solution 
of Eqs. (2.38|)-(2.39) can be expressed in the following form: 



X(t) = 2 V —{[q(t),q(t')]Z(t,t',t )} 



'dt 




y(t)=2r)—{[p(t),q(t')]E(t,t',t )} , 



(2.40) 
(2.41) 

(2.42) 

where 9(t) is the Heavyside step function (defined so that 8(0) = 1/2). Using Eqs. ( [2.36 ) within Eq. ( 2.40 ) and setting 
t = t' , it is possible to observe that 



with 



S(t, t', t ) = 9(t - t Q )6(t - t')6{t' - t ) - e(t - t)6{t' - t)9(t - t'), 



X(t') = [q(t'),Q(t')] = ls(t'-t )[q(t'),p(t')] 
J m 



(2.43) 



Finally, using this results in Eq. ( 2.37| ) with t' — t, one gets the desired result, i.e., that the time derivative of the 
commutator between q(t) and p(t) is equal to zero. Similar arguments can be used to prove the preservation of all 
the other commutation relations between Brownian particle operators. 
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III. HOMODYNE SPECTRUM 



Let us now consider again the dynamics of the optical mode of the cavity with a movable mirror. Using the results 
of the preceding section, the dynamics for t > to of the system can be described by the fol lowi ng set of co upled 
quantum Langevin equations in the interaction picture with respect to hcoob^b (see Eqs. ( [2.1| ), ( |2.3| ), and ( 2.36| )), 



q(t) =p(t)/m 



huj r 



p(t) = -mu 2 s q{t) + Q(t) - -Lp(t) + -jftf(t)b(t) 
b(t) = - (iw c - iljo + y ) b{t) + i^-q(t)b(t) + E - 



(3.1) 
(3.2) 
(3.3) 



where the commutation relations and correlation functions of the input noise bi n (t) are given respectively by Eqs. ( |2.4| ), 
(U) and (U|), while those of Q(t) are given by Eq. fl2.29Q and Eqs. ( |2.32| )-( |2.34| ) considering the limit -> oo. 

In standard interferometric applications, the driving field is very intense. Under this condition the system is 
characterized by a semiclassical steady state with the internal cavity mode in a coherent state \B st ), and a new 
equilibrium position for the mirror, displaced by q s t = TiLo c \B st \ 2 /maj^L with respect to that with no driving field. 
The steady state amplitude is given by the solution of the nonlinear equation 



Be* — 



E 



-if + lUc 



ILOq 



(3.4) 



which is obtained by ta king the expectation values of Eqs. (3.1 )-(3.3) , factorizing them and setting all the time 
derivatives to zero. Eq. (3.4) shows a bistable behaviour which has been experimentally observed in [ B6fl . 

U n der t hese semiclassical conditions, the dynamics is well described by linearizing the quantum Langevin equations 
( 3. 1 )- ( 3.3 ) around the steady state. If we now rename with q(t) and b(t) the operators describing the quantum 
fluctuations around the classical steady state, one gets 



q(t) = p(t)/m 



1 hujAt (b(t) + bHt))+Q(t) 



p{t) = -mu%q{t) - -Lp(t) + 
m 



L 

c + iA) bit) + i^y^q(t) + VTM) , 



V2 -y— - L 

where we have chosen the phase of the cavity mode field so that B st is real and 



(3.5) 
(3.6) 

(3.7) 



A = ui c - wo - 



fiwl 



\B st \ 2 



(3.8) 



is the cavity mode detuning. We shall consider from now on A = 0, which corresponds to the most common 
experimental situation, and which can always be achieved by appropriately adjusting the driving field frequency luq. 
In this case the dynamics becomes simpler, and it is easy to see that only the phase quadrature Y(t) — i (b^(t) — b(t)) 
is affected by the mirror position fluctuations q(t), while the amplitude field quadrature X(t) = b(t) + b'(t) is not. In 
particular, in the limit of a sufficiently large cavity mode bandwidth 7 C 3> 77/774,075 (which is usually satisfied), the 
dynamics of the phase quadrature adiabatically follows that of the mirror position, that is 



Y(t) = ^q(t) 



noise terms. 



(3.9) 



Therefore a phase noise measurement, as for example the homodyne measurement of the field quadrature Y(t), gives 
a direct information on the mirror Brownian motion. In particular, the interesting measurable quantity is the output 
power density spectrum of the phase quadrature 



Sy(w) 



2tt 



dTe luJT (Y out (t)Y out (t + T)) 



dul ' e ^+^')t { Y out (cj')Y out (cj)) \ , 



(3.10) 
(3.11) 
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where {•••}* denotes the time average over t, Y out (t) — i b out (t) — b out (t) is the output phase quadrature, Y out (ui) 
f dte tu>t Y out (t) is its Fourier transform, and 

b ut(t)+b m (t) = V7c6(*) (3.12) 



is the usual input-output relation j27|. If we now take the Fourier transform of Eqs. ( |3.5| )-(3.7), one gets the following 
expression for the Fourier transform of the correlation function of the cavity mode annihilation operator 



&(«) 



1 



D(lo) 



u c B s t ( huj c B sty /%b in (uj) + b\ n {uo) 



mL 



L 



iio - 7 c /2 



Q(lo) - (lo% - to 2 - irjw/m) ^/i^b in (uj) 



with 



D( 



(lo) — (iio — —J (uJg — lo 2 — ir/Lo/m) 



(3.13) 



(3.14) 



For the evaluation of the phase noise spectrum Sy(lo), one has to use Eq. (3.13) and then consider the spectrum of 
t he va ri ous n oise t erms , which are easily derived from the Fourier transform of the correlation functions (2.5), ( ^.(| ). 
( p2| ), fl2~33l ) and (f04|): 



(Q(lo)Q(lo')) = 2nr]huj 



1 + coth 



\2kT 



5(UJ + Lo') 



{b\ n {L0)b m {u')) = 

(b in (uj)bt n (uj'))=2Trd{u + u') , 

where we have assumed again the infinite cutoff limit Q 
operator Q{t). One finally obtains 



(3.15) 

(3.16) 
(3.17) 



oo in the evaluation of the spectrum of the thermal noise 



S Y {w) = 1 + 4 



^ tlLU 2 ^ c \B s 



iL 2 



1 



(( 7c /2)2 +w 2)|DH|2 



u 2 c Vlc\B s 
m 2 L 2 



\D(u)\ 



- Tilj coth 



/ Tiuj 



\2kT 



(3.18) 



This is the phase noise spectrum associated with the homodyne measurement of the phase quadrature Y(t), and 
the only assumptions made in its derivation are the linearization around the semiclassical steady state and the time 
coarse-grained description fl — > oo. Its temperature dependence is instead exact and therefore Eq. ( 3.18 ) is valid even 
at very low temperatures, differently from the spectra obtained with the approaches based on the master equation, 
as in Ref. [jT2"[ , which cannot be applied in the low temperature limit. 

Notice that SV (o>) is an even function of lo, as it must be due to stationarity and the commutation rules of output 
fields ]l2[] . In fact, the non-sensical term of the spectrum found in Ref. (l2) in the case of the standard quantum 
Langevin description is due to the inappropriate form of the correlation function of the Langevin nois e dic tated by 
the SBMME, and it is absent when the correct spectrum of the quantum Brownian noise term of Eq. (3.15) is used. 
This spectrum is plotted in Fig. 2 (see the figure caption for parameter values), where also the three contributions to 
the noise spectrum are explicitely shown. The full line refers to the total homodyne spectrum, while the dashed-dotted 
line describes the shot noise, which is frequency-independent (actually Sy(lu) has been defined so to be normalized 
just to the shot noise level). The dashed line describes the second term in Eq. (3.18), which is the one associated with 
the radiation pressure; finally the dotted line describes the last term which is just the thermal noise contribution. 

The ho mody ne spectrum derived in Ref. jl2| with the adoption of the Diosi master equ ation of Ref. |l4|] coincides 
with Eq. ( |3.18 ) except for a different thermal noise term, which is obtained from that of Eq. (3.18) with the replacement 



huj coth 



Tvuj 
2kf 



2kT + 



[to" 



2 1^,2 

rf jm 



6kT 



(3.19) 



However, despite this formal difference, the two predictions become practically indistinguishable if typical experimental 
parameters are considered. I n fac t, the prediction of Ref. [jl2| coincides with the high temperature expansion (at 
first order in Tioj/kT) of Eq. ( 3.18| ) except for the additional factor firj 2 /6m 2 kT '. However, in typical experiments, 
mechanical oscillators with a very good quality facto r are always used, so that the term h 2 r) 2 /6m 2 kT will be in practice 
always negligible with respect to H 2 lo 2 /6kT in Eq. ( 3.19| ). This means that an appreciable discrepancy between the 
two expressions of the thermal noise term manifests itself only when kT < Hlo, which means prohibitively small 
temperatures, or alternatively, very large frequencies (larger than 1 THz at liquid He temperatures). Moroever at 
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these high frequencies the thermal noise contribution is completely blurred b y the shot noise term and therefore we 
can conclude that with present tecnology the phase noise spectrum of Eq. (3.18 ) and that evaluated in Ref. 
cannot be experimentally distinguished. Nonetheless, the result of Eq. (3.18) is important because it s hows that 
the standard quantu m La n gevin equations (supplemented with the appropriate commutation relations ( 2.29 ) and 
correlation functions ( 2.32 )-( 2.34 ) of the random Langevin force) do give an adequate description of quantum Brownian 
motion, which is even more general than that associated with the master equation, which is not valid at very low 
temperatures. 



IV. CONCLUSIONS 



We have considered in this paper the dynamics of a cavity mode with a movable mirror, which is often used for 
the interfcrometric detection of very weak forces. We have focused in particular on the description of the quantum 
Brownian motion of the mirror, which is responsible for the thermal noise term in the measured phase noise spe ct rum 
of the light reflected from the cavity. We have shown that the standard quantum Langevin equations (fO|)-([2~9|) 
provide an adequate and consistent description of quantum Brownian motion. We have derived the quantum Langevin 
equations directly from the independent oscillator model (providing the commonly used description for the oscillator 
reservoir, see [jl3|,|l5],^3 24 ), and we have seen that they provide a quite general description of quantum Brownian 
motion, valid at any temperatures. This is instead not true for master equation-based approaches, which cannot be 
applied in the low-temperature limit p2| | . The inadequacies found in the quantum Langevin approach are to be traced 
back to the fact that the quantum Langevin force appearing in it is different from the standard input noise terms of 
the input-output formalism p5| , since it is characterized by a different commutation relation (see Eq. (2.29)) which 
does not coincide with a Dirac delta in any limit. 



APPENDIX A: 



In order to justify the presence of the sign function S(t) on Eq. ( 2.36 ), let us consider the step function 



1 if r > 

9(t) = { 1/2 if r = (Al) 
if r < , 



and the formal identity: 

l(t,t ) = I dt'5(t-t')p(t') 



to 
oo 



dt'S(t - t')p(t')(9(t' - t ) - 6(t' - t)) (A2) 

which holds for every t and to. Now using the formal properties of S(t) on Eq. ([A2]), one can verify that for t > to and 
for t < t , T(t,t ) = p(t)/2 and l(t,t ) — —p(t)/2 respectively, while, of course, l(to,to) = 0. This can be written 
using the sign function S(t) — 8(t) — 0(—t) in the following way 

l(t,t )=S(t-to)p(t)/2. (A3) 
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LOCAL OSCILLATOR 



LASER 




FIG. 1. Schematical description of the system. The cavity mode is driven by the laser which, thanks to the beam splitter 
BS, provides also the local oscillator for the homodyne measurement of the light reflected by the cavity. 
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FIG. 2. P hase noise spectrum of Eq. ( 3.1S ) (full line). The dotted-dashed line refers to t he sh ot noise spectrum (first term 
of Eq. (3.15)); the dashed line refers t o the radiation pressure term (second term of Eq. (3.18)), and the dotted line to the 
thermal noise term (third term of Eq. (3.18)). Parameter values are los = 1.3 ■ 10 5 Hz, rt/m = 3 ■ 1(T 2 Hz, lo c = 1.8 • 10 15 Hz, 
7 C = 4.7 ■ 10 5 Hz, m = 10" 5 Kg, L = lCT^m, T = 4.2 K, P = 10" 5 W. 
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